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Abstract 

We reconsider the structure-based route to coarse graining in which the coarse-grained model is 
defined in such a way to reproduce some distributions functions of the original system as accurately 
as possible. We elucidate why standard expressions for pressure and chemical potential applied to 
this family of coarse-grained models with density-dependent interactions fail to reproduce pressure 
and chemical potential of the underlying original system, under the only assumption that the 
interactions in the underlying model be state independent. Moreover, we show that the state- 
dependent potentials depend on the ensemble in which they have been derived. Therefore, care 
must be used in applying canonical state-dependent potentials to predict phase lines, which is 
typically performed in other ensembles. 
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I. INTRODUCTION 



In condensed matter physics, chemistry, and material science state-dependent interactions 
arise in many different contexts. They are usually obtained from coarse-graining procedures 
in which a subset of the degrees of freedom is integrated out (see Refs. ll|-|4] for an overview 
of methods and applications). Indeed, any coarse graining of the original (be it classical or 
quantum) system induces many-body interactions among the remaining degrees of freedom. 
The idea is then to replace these complex interactions with state-dependent potentials that 
are usually taken to be pairwise additive for computational efficiency (if necessary, many- 
body forces can also be introduced, as in Ref. Q) and that are therefore much more tractable 
from a theoretical and/or numerical point of view. Different criteria have been used to select 
the optimal set of state-dependent pair potentials. In the structural route the model with 
state-dependent interactions is required to reproduce some distribution functions associated 
with the coarse-grained (CG) degrees of freedom.- - — It is also possible to define the CG 
model by matching the forces^— acting on the CG sites computed in the original, underlying 
system, or by requiring the CG model to reproduce solvation free energies. Also state- 
dependent potentials suitable to treat inhomogeneous systems have been proposed.— 1 ^ We 
should further mention mixed coarse-graining strategies that try to match simultaneously 
the pair distribution function and some other thermodynamic property, for instance, by 
constraining the virial pressure to be equal to the pressure of the microscopic model.- - - 1 ^ 
These approaches, although in principle incorrect since the potential is uniquely defined 
by the pair distribution function g(r) and the density according to Henderson's theorem,— 
may still be of value in practical numerical calculations. Indeed, g(r) is only known with 
statistical errors and is little sensitive to the tail of the potential: visibly different potentials 
may produce structures with essentially identical pair distribution functions.— Therefore, 
the large-distance behavior of the CG interactions is determined with a large uncertainty, 
which might leave some flexibility to implement an additional constraint. We will not discuss 
these approaches any further, focusing on the conceptual problems of the approach rather 
than on the difficulties of practical numerical implementations. 

In this paper we discuss the structural approach, which dates back to the early days 
of liquid state theory,- - — considering state-dependent pair interactions for a very general 
class of classical systems. The microscopic, underlying model to which the coarse-graining 



2 



procedure is applied consists of polyatomic molecules. No hypothesis is made on the nature 
of the interactions among the atoms: we only require them to be state independent, but, 
otherwise, they are arbitrary and can, in principle, include any type of many-body terms. 
Our discussion therefore applies both to monoatomic systems with three- and higher-body 
interactions (for instance, to noble gases whose accurate study requires the introduction of 
the three-body Axilrod- Teller potential) and to soft-matter systems, such as polymers, 
in which the complexity lies in the number of atoms involved rather than in the many- 
body nature of the atom-atom interactions.— 1 ^ 1 ^ Within the CG approach we replace the 
microscopic system with a CG system of monoatomic molecules that interact by means 
of state-dependent pair potentials. The latter are fixed by a structural requirement, the 
equality of a specific pair distribution function. If the underlying system is formed by 
monoatomic molecules, we consider the usual radial pair distribution function. In the case 
of polyatomic systems, each molecule is replaced by a CG monoatomic molecule located 
at some point X. In many instances, X corresponds to the center of mass of the original 
molecule, but other choices are possible: for instance, in CG models for star polymers the 
point X usually coincides with the center of the star.— ~— For our purposes, we do not need 
to specify how X is chosen. We only require the coordinates of X to be a weighted 
average of the positions of the atoms belonging to the molecule. Once X is chosen, we can 
consider the X-X pair distribution function 



where p is the density and N the number of molecules. The state-dependent potential is fixed 
by the requirement that the pair distribution function gco{ r ) m the CG model is equal to g(r) 
in the underlying system. From a practical point of view, the determination of the potential is 
not an easy task and several method have been devised, like the iterative Boltzmann inversion 
method^ and the inverse Monte Carlo method.— 1 ^ Variational approaches have also been 
discussed, optimizing the coefficients of suitable parametrizations of the CG potentials.— 1 ^ 
Note that the structural approach we have described here is not the only one that is used 
in practical applications. We should also mention the force-matching approach (often called 
multiscale coarse-graining method^), in which the state-dependent potential is determined 
by requiring the CG system to match the atomistic force on the CG atoms as accurately 
as possible. Also this method has a structural interpretation: the matching condition is 




(1) 
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equivalent 16 to require the CG force to satisfy the appropriate Yvon-Born-Green equation 32 
that relates the pair and the three-body correlation function. We mention that general 
theoretical formalisms within which different CG approaches can be derived are presented 



in Refs. 
see Ref . 



33 
35 



34 



For a summary of methods and a detailed comparison in a specific example, 



Structurally derived state-dependent potentials have been mostly discussed in the context 



of the canonical ensemble, see Refs. 



36-38 



and references therein. In this case, the potentials 
become density and temperature dependent. For this reason, some thermodynamic identities 
are no longer satisfied: in particular, the compressibility route and the virial route to the 
pressure are no longer equivalent for the CG system. This gave rise to several investigations 
concerning the thermodynamic consistency of models with state-dependent potentials.—* 3 ^ 
In this paper we again investigate the conceptual issues which arise when considering state- 
dependent interactions. The basic question we wish to address is whether the knowledge of 
the state-dependent potentials for a given value of the density p allows one to determine the 
exact value of thermodynamic quantities, such as chemical potential and pressure, for the 
underlying system at the same value of p using the standard tools of statistical mechanics. 
Here, we mainly extend the discussion of Ref. y. We show in full generality that the virial 
pressure or the generalized expression of Ascarelli and Harrison 3 - do not provide the correct 
underlying-model value. The same result holds for the chemical potential: Widom's method 
does not provide the correct answer. Note that, if we assume that the density-dependent 
(DD) potentials are known for all densities < p < p max and no phase transitions occur in 
such density interval, then, via the compressibility route, we can compute the exact equation 
of state of the underlying model for p < p max , hence derive the exact pressure, chemical 
potential, and free energy. However, in this case DD potentials play little role. Indeed, 
their computation requires a study of the microscopic model in the whole density interval. 
But this study would also provide the equation of state directly, without the necessity of 
introducing any sort of CG model. 

It is interesting to extend the analysis of the thermodynamic properties of state-dependent 
potentials to other ensembles. For instance, for phase-coexistence studies it is more natural 
to consider the grand-canonical ensemble, while the isobaric ensemble may be the best suited 
to interpret experimental data, since pressure is fixed in experiments. In these ensembles one 
would consider fugacity- and pressure- dependent pair potentials, respectively. Again, the 



4 



question we wish to address is how to obtain thermodynamic predictions in these different 
ensembles by using state-dependent interactions. 

The paper is organized as follows. In Sec. [IT] we discuss density-dependent potentials 
and the consistency of several commonly used methods to determine pressure and chemical 
potential. We show in full generality that none of them is exact: none of the thermodynamic 
quantities extracted via a CG approach is more than an approximate estimate. In Sec. II III 
we study the accuracy of the different approximations: first, in Sec. IIII Al we present a 
general, model-independent discussion in the low-density limit, then, in Sec. IIII Bl we present 
a specific calculation for linear polymers under good-solvent conditions. The analysis of 
Sec. HUis extended to the grand-canonical ensemble in Sec. HVl while in Sec. |V]we discuss the 
accuracy with which the thermodynamic behavior is reproduced by using fugacity-dependent 
potentials. Finally, in Sec. [VTJwe present our conclusions. 



II. DENSITY-DEPENDENT POTENTIALS IN THE CANONICAL ENSEMBLE 

In the structural approach we consider, the basic quantity of interest is the X-X pair 
distribution function g(r; p) defined in Eq. (TjQ) for the underlying model. Knowledge of 
g(r; p) allows us to obtain the thermodynamics of the system by using the compressibility 
relation,™ which holds for any choice of the point X: 

Y^ = l + pj d 3 r[g(v;p)-l], (2) 

where 

K{p) = (3) 

is related to the isothermal compressibility xt = P/(pK) an d (3 = 1/ksT. In the absence of 
phase transitions in the density interval [0, p], other thermodynamic quantities like pressure, 
chemical potential, and Helmholtz free energy at density p can then be obtained from K(p) 

as 

0P{p) = [ P daK(a), (4) 
Jo 

f3p(p) = ln^+ f P ^(K(a)-l), (5) 
-]y = 0f(p) = /W) — > (6) 



where q\ = Z\/V and Z\ is the partition function of a single, isolated molecule. For conve- 
nience, here and in the following we set the de Broglie thermal length equal to one. 

In the CG approach we discuss here, one replaces each molecule with a monoatomic one. 
The CG molecules interact by means of the DD potential Vdcg( t i P)i which is defined such 
as to reproduce the finite-density pair distribution function g(r; p). If we consider the density 
appearing in the pair potential as a parameter, we can define a standard thermodynamic 



system (this is the passive approach of Ref. l36l ). Therefore, we define the Helmholtz free 
energy as 

F CG (N,V,T;p p ) = Nf CG (p;p p ) (7) 
= — A^Tln J dr x ... dr N e"^ E « ^coin-r^) ^ 

where p = N/V. The free energy (J7J) depends on two densities: p = N/V is the usual 
quantity, while p p is the density parametrizing the pair potential. Using Eq. ([7]) we can 
define the pressure and the chemical potential for the CG model at fixed p p : 



UV ' N,T,p p 



,caM={^) vTp . (9) 

If gcG{ r ', Pi Pp) is the pair distribution function in the CG model with free energy ((7j), the 
compressibility relation (jSJ) allows us to compute K GG (p, p p ). Then, the standard thermo- 
dynamic relations give us 

Ppcg(p,Pp) = lnp+ / — (K CG (cr,p p ) -1), (10) 

Jo ° 

PPcg(p,Pp) = / da K CG (a,p p ), (11) 
Jo 

Pfaofa Pp) = PPcg(p, P P ) ~ (3PcG ^ Pp \ (12) 

These expressions are the analog of those appearing in Eqs. (jl]), (j3j), and (jSJ), with q± = 1, 
since the CG model we consider consists of monoatomic molecules. Since p p is a fixed 
parameter, the pressure P G g{p, Pp) an d the chemical potential Pcg{p, Pp) can be determined 
by using standard methods. The virial expression holds for the pressure, hence we have^ 

I3Pcg(p,P p ) = l3Pviv(p,P P ) = P 

-^f «^ g(w ,A, (13, 



Analogously, the chemical potential can be obtained by using Widom's insertion method ■ 



/3pcg(p,P p ) = lup- 



in 



i j d 3 r N+1 ( e -PU N+1 (.r N+1 ;Pp)) 



N,V 



(14) 



where (-)n,v is the canonical ensemble average over N molecules in a volume V, p — N/V, 
and Un+i ( r N+i ] Pp) is the insertion energy of an additional molecule computed by using 
V D , CG (r;p p ). 

The DD pair potential at density p p = p is determined by requiring the pair distribution 
function gccij] Pi p) to be equal to the XX pair distribution function g(r; p) in the under- 
lying microscopic model: g GG (r; p, P) — g( r 'iP)- Because of the compressibility relation (J2J), 
this implies the equality of the compressibilities, i.e. 

K{p) = K CG {p,p). (15) 

This equality does not extend, however, to the other thermodynamic quantities: if we set 
p p = p in pcG or PcGi we do not obtain the correct chemical potential and pressure at 
density p in the microscopic model. Indeed, barring unexpected concidences, we expect in 
general 

(3p exc (p)= [ P — (K CG (a,a)-l) 
Jo <y 

^(3p^ G (p,p)= [ P ^(K CG (a,p)-l), (16) 

(3P(p)= f ' daK CG {a,a) 
Jo 

± (3P CG (p,p) = [ P daK CG (a,p), (17) 

where we used Eq. (TT5J to replace K(p) with K C g(p, p) in the left-hand side integrals, which 
give chemical potential and pressure in the underlying model. Hence, if we use the virial 
route for the pressure or Widom's insertion method for the chemical potential in the CG 
model, we do not obtain the exact result, but only approximations to the correct pressure 
and chemical potential: even with the use of DD potentials, we cannot obtain the correct 
thermodynamics at density p from the study of the CG model at that density. Of course, 
the pressure and chemical potential of the underlying system can be correctly obtained by 
using K(p) = K C g{Pi p) an d relations (j4j) and (JSJ). But this approach essentially uses the 
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equation of state of the microscopic model and does not make use of the CG model with 
DD potentials. 

In the context of DD potentials, Ascarelli and Harrison^ proposed a different formula 
for the pressure. The idea is to start from Eq. (j7j), set p p = p = N/V, and then take the 
derivative with respect to V (this is the active approach of Ref. 1361 ) . This gives 



Pah(p) = P 



2 dfcc(p, p) 



dp 

Pca(p,p) + P 2 
Pvh(p,p) + 



dfcG(p,P P ) 
dp p 



p p =p 



(18) 



2vrp 3 



r dr ■ 



dp 



where P v i r (p,p) is explicitly given in Eq. ( fl3l) . Simple algebra allows us to rewrite 



PPah(p) 



da 



K C g{o, p) + 



P(P ~ °) dK CG (a, p) 



(19) 



a dp 

which shows that Pah(p) differs in general from the pressure of the underlying model. 
The Ascarelli-Harrison prescription can also be applied to the chemical potential, defining 

9fcG(p,p) 



Vah(p) = fcc(p,p)+P 

= Vcg{p,p)+P 

= Vcg{p,p) + 
2vrp 2 



dp 

dfcoipiPp) 



dp P 



p P =p 



(20) 



r dr 



dVi 



D,CG 



>;p) 



dp 



g CG (r;p,p)- 



Again, it is easy to verify that such an expression does not give the correct result. Indeed, 



^ah(p) = 
which differs from f3p exc (p). 



p da 
a 



K CG (cr,p) -l + (p-a 



dK CG (a, p) 
dp 



(21) 



III. COMPARING THE DIFFERENT DEFINITIONS OF THE PRESSURE AND 
CHEMICAL POTENTIAL IN THE CANONICAL ENSEMBLE 

As we have seen in the previous section, the DD quantities Pcg(p, p) an d pcg(p, p) 
provide only approximations to the correct result. In this section we wish to compare them 



S 



with the exact result, determining quantitatively the size of the discrepancy. In the context 
of soft-matter systems, it is also common to use CG models based on zero-density potentials, 
i.e. potentials obtained by matching the pair distribution function in the limit p — > 0. Here 
we will compute the corresponding pressure Pcg(Pi 0) and chemical potential Pcg(Pi 0) and 
we will compare them with the corresponding DD quantities, in order to understand which 
method provides the best approximation. We shall first consider the low-density limit, 
which can be analyzed in a model-independent way, and then we shall apply all formulae to 
a specific soft-matter example, polymers in the semidilute regime. 

A. Low-density behavior 

Let us consider the low-density limit. For p — » the XX pair distribution function can 
be expanded as 



and A(r) = g\(r) — ^i(r). The quantity A(r) encodes the contributions of the three-body 
interactions: for monoatomic systems interacting by means of pairwise additive interactions 
we have A(r) = and <7i(r) = <?i(r). Using the compressibility relation (jSJ) and Eq. (1221) . we 
can compute K(p) in the underlying microscopic model: 



#(r; p) = 00 (r) + P9i( r ) + 0(p 2 ). 



(22) 



We define h (r) = g (r) — 1, 




(23) 



K(p) = l-pI -p 2 (I 1 + I 2 ) + O(p 3 ), 



(24) 



with 




(25) 



(26) 



(27) 



Eqs. (jU) and (j^D give 



f3p e * c (p) 



PP(P) 



P-^ 2 /o-^P 3 (/i + / 2 ) + 0( P 4 ), 
-I p-^(I 1 + I 2 ) + O(p 3 ). 



(28) 



(29) 
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Let us now consider the CG model. First, we need to compute the DD potential Vd,cg( t ', Pp)- 
Since we are interested in the low-density limit, we expand it as 

Vb,cc(r; Pp) = V c G (r) + p p V WG (r) + 0{p 2 p ). (30) 

At leading order in the density, we have 

g CG (r;p,p P ) = C -*W*> + 0(p )Pp ). (31) 

By requiring the equality of the pair distribution functions, i.e. g(r; p) = gccij] Pi p)i we 
obtain 

0V OCG (r) = -ln^o(r). (32) 

If we include the first density correction we have 

g CG (r;p,p p )=e-^ c ^[l+ (33) 
p J d 3 sf(s;p p )f(s-r;p p ) + 0(p 2 ), 

where 

/(s;p p )= e -^^( s ^-l. 
Expanding in p p , using Eq. (152"]) and h (r) = g (r) — 1, we obtain 

gcc(r; p, pp) = go(r) [1 - p P l3V 1CG (r)] 

+pg (r) J d 3 sh (s)h (r-s) + O(p 2 ,pp p ,p 2 p ). (34) 

By requiring g(r; p) = gcG^] P, p), we obtain the next correction: 

PVwaW = -^4- (35) 

Note that these definitions imply 

J d 3 v /(r; p p ) = I + p p I 2 + 0(p 2 p ), (36) 

for p p — > 0. We can now compute K(p,p p ). In the limit p — > 0, using the standard 
expressions valid for a monoatomic system, 32 we obtain 

K CG (P, P P ) = 1 - P J d 3 r /(r; p p ) (37) 
-p 2 J d 3 rd 3 s /(r; p p )f{s; p p )f{r - s; p p ) + 0(p 3 ). 

10 



Expanding in p p and using Eq. (1361) . we obtain 

K C g(p, P p ) = 1 - ph ~ P 2 h ~ PPph + 0(p 3 , p 2 p p , pp 2 p ). (38) 

For p p = p this expression coincides with Eq. (T2"4"|) . as expected. Using Eqs. ffTUj) . (fTT|) . and 
(|T2|) we obtain 

f3P CG (p, pp) = p- \p 2 h ~ \p*h ~ \p 2 p v K (39) 

PPcg(PiPp) = ~Ph ~ \p 2 h - PP P h, (40) 

Pfccip, P P ) = hip - 1 - ^pl - ^p 2 h - ^PPph, (41) 

disregarding higher-order terms in p and p p . From the expression of the free energy, we 
obtain the Ascarelli-Harrison expressions for the pressure and chemical potential: 

0PaM = p/J^± = p - \p% - \p% - p%, (42) 

PpTh(p) = Pfcc(p,p)+Pp dfcG ^ p p) 

= -Ph ~ \p 2 h - \p 2 h- (43) 

We can thus compare these expressions with the exact ones valid for the underlying model. 
For the pressure we have 

(3P CG (p,p)-(3P(p) = -lp 2 I 2 , 

o 

(3P CG (p,0)-(3P(p) = ^p 2 I 2 , 

(3P AH (p)-PP(p) = -lp 2 h- (44) 
From these results we see that the pressure computed by using the DD potentials provides 



the best approximation to P{p)- The AH formula is significantly worse — deviations are 



37 



our 



38 



times larger than those for the virial pressure Pcg(p, p) — as already noted in Refs 
By using the zero-density potentials one obtains a result which is worse by a factor of two 
than the DD result. 

We can perform the same comparison for the chemical potential. We find 

PPcciPi P) - PP EXC (P) = -\p 2 h, (45) 

/3^(p,0)-/3/i exc (p) = ip 2 / 2 , (46) 

PpTh(p) - PP SXC (P) = -P 2 h- (47) 
11 



TABLE I: Compressibility factor Z = f3P/ ' p and excess chemical potential for a polymer CG 
system at $ = 0.8 and 1.5. The reference polymer Z(p) and p cxc (p) have been obtained by using 
the equation of state reported in Ref. Q- 

3> Z{p) Z CG (P,0) Z CG (p,p) Z AH (p) (3p™(p) gMg^O) /VcgOM 
0.80 2.35 2.22402(4) 2.44790(4) 2.51(2) 2.56 2.3797(2) 2.8108(3) 
1.50 3.90 3.39847(4) 4.13973(4) 4.28(5) 5.38 4.653(4) 6.0929(4) 

In this case, the chemical potential derived by using the DD and the density-independent 
potentials have the same accuracy. Expression (|20|) is instead worse than both the one 
obtained by using DD potentials and the one obtained by using zero-density potentials. 

Note that all deviations are proportional to I2, which encodes the contributions of the 
(effective) three-body interactions. Hence, consistency is only possible if the three-body 
terms do not contribute to the thermodynamics, hence I2 = 0, a result which is not generally 
true since the free energy of the CG model in a generic CG procedure is never pairwise 
additive. 



B. Polymers in the semidilute regime 

To assess the quality of the different approximations to the pressure and the chemical 
potential, we study the thermodynamic behavior of a system of linear polymers under good- 
solvent conditions. We have extensively studied CG models for this type of systems in 



Refs. 



41 



42 



For two values of the polymer volume fraction $, 

47T 



(48) 



where R g is the zero-density radius of gyration, $ = 0.80 and $ = 1.50, we compute^ 
the center-of-mass pair distribution function g(r;p). Then, by using the method described 
44j, we determine the DD potential Vd i cg( i ';p); checking that it reproduces the 



in Refs. 
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full-monomer value of K. This is a strong check of the accuracy with which the potential 
is determined, since the compressibility is very sensitive to the large-distance behavior of 
the pair potentials.— For the zero-density potential we use the accurate expression reported 



in Ref. 



45 



The Monte Carlo estimates of the compressibility factor Z = (3P/p and of the 



chemical potential obtained by using the virial route and Widom's method, respectively, are 
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reported in Table [fl Completely consistent results are obtained by using the integral-equation 
approach and the hypernetted-chain (HNC) closure. In our previous work— we showed that 
the compressibility factor Z C g(p, 0) obtained by using zero-density potentials underestimates 
the polymer Z(p), which, in the low-density limit, implies I2 < 0. Our present results are 
fully consistent. Differences increase with increasing $, as expected: the relative error is 5% 
and 13% for $ = 0.8 and 1.5, respectively. If one uses the DD potentials the compressibility 
factor is overestimated, by 4% and 6% in the two cases. Therefore, in agreement with 
the low-density analysis, DD potentials provide more accurate estimates of the pressure 
and clearly represent an improvement with respect to zero-density potentials, especially 
for the largest value of $. We also compute the pressure by using the Ascarelli-Harrison 
expression. A precise determination is not easy: indeed, it is difficult to estimate accurately 
the derivative dVD,cG( r ] p)/dp, especially for r/R g > 2. Hence, the additional term that 
appears in the Ascarelli-Harrison expression can only be determined with limited precision, 
which we estimated somewhat subjectively by looking at the variation of the results when 
different parametrizations of the DD potentials are used: we estimate a relative error on the 
correction term of approximately 20-30%. As expected, Zah is significantly worse than the 
virial expression Zcg{p,p)'- the relative error is 7% and 12% for the two values of We 
also considered the combination P'(p) = ^Pcg(p, p) — ^Pah(p), which is correct to order p 2 , 
see Eq. (jS]). We find Z'{p) = 2.43(1), 4.06(2) for the two values of $. These estimates are 
slightly closer to the exact result than Zcg{p,p)- 

Similar conclusions are reached for the chemical potential. If one uses zero-density po- 
tentials, /3/i exc (p) is underestimated by 7% and 14% in the two cases, while, by using DD 
potentials one overestimates the chemical potential by 9% and 13%. As already discussed 
in the low-density regime, zero-density and DD potentials are equally inaccurate. 

This discussion indicates that, for $ < 1, which is the region in which CG monoatomic 
models should work reasonably, there is little advantage in using DD potentials. Accurate 
results can only obtained by using multiblob models.— 4241-^2 For instance, the tetramer 



model (model 4MB-2) of Ref. 



42 



gives Z(p) = 2.3505(4) and f3p exc (p) = 2.57(1) for $ = 0.8, 
in perfect agreement with the full-monomer results. For $ = 1.5, small deviations, of the 
order of 2-3%, are found: Z = 3.7901(2) and (3p exc = 5.27(2). If more accuracy is needed, 
one can increase the number of blobs. If one uses a decamer with 10 blobs,— one obtains 
for $ = 1.5 the estimates Z = 3.929(3) and /3p exc = 5.3(1), which are consistent with the 
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full monomer results. 



IV. FUGACITY-DEPENDENT POTENTIALS IN THE GRAND-CANONICAL 
ENSEMBLE 

The previous discussion applies to the canonical ensemble. However, if one is interested 
in studying phase coexistence, the grand canonical ensemble is the natural one. In order to 
implement the CG procedure, we parametrize all thermodynamic variables in terms of the 
reduced fugacity z = fte^*, where q\ = Zi/V and Z\ is the partition function of a single 
molecule. This guarantees that z ~ p for p — > 0, both in the original and in the CG model. 
As before, the basic quantity in the approach is the XX pair correlation function g(r;z), 
which is now a function of the fugacity z. If p(z) gives the density of the underlying system as 
a function of z, we have obviously g(r; z) = g(r; p(z)). The relation with the thermodynamics 
is provided by the compressibility relation. In the grand canonical ensemble the role of K 
is played by 

which satisfies^ 

H(z) = z J dh[g(j-z)-l}. (50) 

Using the standard thermodynamic relations we obtain (of course, we assume here that there 
is no phase transition in the interval [0, z\) 



P{z) 



z dw n ~ 



1- / — H(w) 
o w 



(51) 



0P{z)= f Z — p(w). (52) 
Jo w 

In the CG model one considers a fugacity-dependent (FD) pair potential V / p i c*G( r ; z )- Again, 
we can consider the fugacity appearing in the potential as a parameter and define 

oo N 
Z 



PQcg{z, Zp) = - In ^2 j^Qn{z p ), (53) 

N=0 

where Qn{ z p) is the partition function of a system of iV monoatomic molecules interacting 
via Vf,cg{ v 'i z p)- The FD potential Vf,cg(?', Zp) is fixed so that the pair distribution function 
in the CG model gcc{ r ] z, z p ) for z p = z is equal to the XX pair distribution function in 
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the underlying model. At fixed z p , we can define the standard thermodynamic quantities: 

P CG (z,z p ) = ~^Q C g(z,z p ), (54) 
P cg(z,z p ) = -- -—— . (55) 



V V dz 

Since Zp IS cL fixed parameter, the density pcg(z,z p ) can be computed as (N) ZiZp /V as 
usual [here is the average with respect to the grand canonical distribution defined by 

Eq. 053p], while Pcg(z, z p ) can be determined by using the grand-canonical generalization 
of the virial expression. If the pair distribution function gcci?', z , z p) is known, we can also 
compute thermodynamic quantities using the compressibility route. First, we determine 

H CG (z, z p ) = z J d 3 r[g CG (r; z, z p ) - 1], (56) 

then Pcg( z , z p) an d (3Pcg(z, z p ) are derived by using Eqs. f[5"Tj) and f[5"2j) at fixed z p . The 

equality of the pair distribution functions implies H C g( z , z ) = H(z), but this result does 

not extend to pc G ( z , z ) and Pcg(z, z). For instance, 

z f z dw TT z f z dw TT , 

? 7=1-/ — Hcg (w,z ^ -rr = 1- / —H CG (w,w). 57 

Pcg(^,2;) Jo w p{z) J w 

It is important to stress that the DD potential Vd,cg{y'-, p) defined in the canonical ensem- 
ble and the FD potential Vf,cg{ v \ z ) are not simply related. If p(z) gives the density as 
a function of z in the underlying model, Vf } cg{ t ) z ) differs from Vd,cg{ t ) p( z ))'- the same 
thermodynamic state in the microscopic model is represented by different state-dependent 
potentials, that depend on the ensemble one chooses. Ensemble equivalence does not hold 
for state-dependent potentials. To understand why the equality does not hold, note that the 
system at fugacity z interacting with pair potential Vf,gg{ y ', z ) has density pcc( z , z ) 7^ p{ z )i 
see Eq. ( [571) . Hence, VF,cG( r ; z) and Vd,cg(?'i p( z )) correspond to two systems that have 
the same pair distribution function but at different densities: hence they cannot be equal. 
The only exception is the zero- density /zero-fugacity case. Since Pcg(0,0) = 0, we have 
Vf^cg^i -2 = 0) = Vd i cg( t ] P — 0). The ensemble inequivalence, which is completely general 
and extends to other ensembles, say the isobaric one, is clearly related to the thermodynamic 
inconsistency of state-dependent potentials, which are not able to reproduce simultaneously 
the correct value of density, chemical potential, and pressure. It is also completely consistent 
with Henderson's theorem,— which states that, if two systems interacting by pairwise po- 
tentials have the same pair distribution function at the same density, then the pair potential 
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is unique. Indeed, in the different ensembles, the pair distribution function is the same, but 
the density is not. 

V. COMPARING THE ESTIMATES OF THE DENSITY AND PRESSURE IN 
THE GRAND CANONICAL ENSEMBLE 

We shall now repeat the analysis of Sec. IHIl in the grand-canonical ensemble, comparing 
the exact p(z) and P(z) with the zero-density approximations pcg(z,0) an d Pcg{z,0) an d 
the approximations pcg{z,z) an d Pcg( z , z ) obtained by using FD potentials. We will first 
consider the low-fugacity limit and then we will give results for a specific example, polymers 
under good-solvent conditions in the semidilute regime. 

A. Low-fugacity limit 

Since g(r; z) = g(r; p(z)) and we have z = p + 0(p 2 ) for p — > 0, we can reexpress Eq. (|22|) 

as 



g{r; z) = g (r) + z[£i(r) + A(r)] + 0{z 2 ). 



(58) 



Using Eq. (150)) we obtain 



H{z) = zlo + z\ll + h + I 2 ) + 0(z 3 ), 



(59) 



where Iq, I\, and I2 are defined in Sec. [TT1 Eqs. (T5T1) and (1521) give 




(61) 



(60) 



Analogously, in the CG model we obtain 



H C g( z , Zp) = z h + z2 (Iq + h) + zz P h + 0(z 3 , z 2 z p , zz 2 ), 



(62) 



which satisfies Hcg( z , z ) = H(z). From this expression we obtain 



Pcg( z , Zp) = z + z 2 I + -z 3 (3Iq + h) + z 2 z p I 2 , 

11 1 

(3P cg ( z , z p) = z + -z 2 I + -z 3 (3I 2 + h) + -z 2 z p I 2 , 



(63) 



(64) 
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disregarding higher-order terms in z and z p . Comparing the CG results with the exact ones 
we obtain for the density 

Pcg(z,z) - p(z) = ^hz 3 , (65) 
p CG (z,0)-p(z) = -\hz\ (66) 

In this case the zero-fugacity and the FD CG models provide approximations that differ 
by the same amount from the correct, underlying- model result. If instead we consider the 
pressure we obtain 

PPc G {z,z)-PP(z) = l -I 2 z\ (67) 

f3P CG (z,0) - (3P(z) = -\hz\ (68) 

o 

It is somewhat surprising that in this case zero-fugacity potentials are more accurate than 
those that depend on the fugacity. We can also compute the compressibility factor Z = 
(3P/p, obtaining 

Z CG (z,z)-Z(z) = -\hz\ (69) 

o 

Z CG (z,0) - Z(z) = h 2 z 2 . (70) 

As in the case of the canonical ensemble, the compressibility factor is instead more accurate 
for state-dependent potentials. 

B. Polymers in the semidilute regime 

We will now repeat in the grand canonical ensemble the computations we have reported 
in Sec. IIIIBI We consider the two values of z, zB? g = 2.47 and 77.7 that correspond to 
$ = 0.8 and $ = 1.5, the values of the polymer volume fraction considered in Sec. IIIIBI To 
compute the FD potentials — as we already noted they differ from the DD potentials — we 



use the same technique employed in Refs. llllJ44j . Since we expect the HNC approximation to 



be quite accurate^ (we verify it below), we relate z and p by using the HNC expressio 



51,52 



n 



pexp 



| J d 3 r (h(r) 2 - h(r)c(r) - 2c(r)) 



(71) 
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where h(r) = g(r; z) — 1 and the direct correlation function c(r) is given by the Ornstein- 
Zernike relation^ 

h(r) = c(r) + p J d 3 s c(s)/i(r - s). (72) 

Given z and g(r;2:), we can solve Eqs. (l7Tj) and f!72|) to obtain p and c(r). The potential 
follows from the HNC closure relation: 

/3Vp jCG (r; *) = h(r) - c(r) - In (/(r; (73) 

In this formalism it is easy to understand why the FD potential at, say, zR? g = 2.47 is 
different from the DD potential at $ = 0.8, although zR 3 = 2.47 and $ = 0.8 correspond 
to the same thermodynamic state of the underlying model. Indeed, Eqs. (ITTj) and f J72|) 
cannot be simultaneously satisfied by fixing the chemical potential and the density to the 
full-monomer values. In the DD case, one fixes p, obtaining a fugacity that differs from the 
full-monomer value, while here we fix z, obtaining a different density. The HNC procedure 
can also be applied to other ensembles. For instance, one could obtain the potential in the 
isobaric ensemble, by simultaneously solving Eqs. f J72|) . f J73|) . and f|T3|) . fixing the pressure 
to its full monomer value. Of course, both p and z would differ from the corresponding 
full-monomer values. The state dependent potentials are reported in Fig. [TJ for the two 
different thermodynamic points of the polymer system we have considered. We also report 
the potential in the isobaric ensemble, obtained by fixing the pressure.— Differences among 
the potentials are not large on the scale of the figure, but still not negligible. At full 
overlap r = 0, we have^ fiVcciT = 0) = 1.775(5) for the zero-density potential (note that 
this potential is independent of the ensemble). At the state point with $ = 0.8 we have 
l3V FtCG (r = 0) » 1.902 for the FD potential and 0V DjC G(r = 0) » 1.955 for the DD potential. 
Hence, in the grand canonical ensemble many-body interactions increase repulsion by 8% 
with respect to the zero-density case. In the canonical ensemble the effect is 40% larger: 
repulsion increases by 11%. At the larger density /fugacity, we have PVp^G^ = 0) ~ 1.955 
and PVd.cg{ t = 0) ~ 2.030 for FD and DD potentials, respectively: repulsion increases by 
11% and 16% in the two cases. Note that the observed differences between the DD and FD 
potentials are of the same order of the difference between the zero-density potential and the 
state-dependent one, i.e. of the order of the contribution of the many-body interactions, 
whose inclusion is the main motivation for considering state-dependent potentials. 

The HNC procedure gives us estimates of $ and of the virial pressure. To check these 
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FIG. 1: We report the DD potential (DD), the fugacity dependent (FD) potential, and the pressure 
dependent (PD) potential as a function of r/R g . They are obtained by fixing: (left) $ = 0.80, 
zR 3 g = 2.47, /3PR 3 = 0.45 in the three cases, respectively; (right) $ = 1.50, zR 3 = 77.7, (3PR 3 g = 
1.40, respectively. 

TABLE II: Compressibility factor Z = /3P/p and polymer volume fraction <3? = AirpR 3 /?> for a 
polymer CG system at two values of the fugacity. The reference polymer Z(z) and &(z) have been 
obtain by using the equation of state reported in Ref. |46|. 

zR 3 g Z(z) Z CG (z,0) Z CG {z,z) <5>cg{z,0) $cg(z,z) 

2.47 2.35 2.2920(1) 2.3338(1) 0.80 0.84104(4) 0.77367(4) 
77.7 3.90 3.7132(2) 3.7819(2) 1.50 1.68500(4) 1.47603(4) 

results, we perform grand canonical simulations of the CG model. We measure the density 
as (N)/V and the pressure by using the standard virial expression. The results reported 
in Table [XT] are in good agreement with those (not reported) obtained by using the HNC 
approximation, confirming the good accuracy of the procedure. 

Let us now compare the results with those of the original, full monomer model. As for 
the compressibility factor, the FD model gives essentially the correct result for the lowest 
value of z (the difference is less than 1%), while it slightly underestimates (by 3%) Z for the 
largest value of z. In both cases, results obtained by using FD potentials are more accurate 
than those obtained by using the corresponding zero-density quantity. The same conclusions 
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are reached for $, especially for the largest value of z. In the analysis of the low-density 
behavior, we observed that FD potentials are worse than zero-density potentials for what 
concerns the pressure. Hence, we also consider 

P = f3PRl = -rZ®. (74) 
y 47r 

At the smallest value of z we obtain P fa 0.46020(3) by using the zero-density potential and 
P = 0.43105(3) by using the FD potential, to be compared with the full-monomer result 
P ~ 0.45. The result obtained by using the zero-density potential is closer to the exact, 
full-monomer result, as expected on the basis of the low-density analysis. At zP? g = 77.7 we 
have instead P ~ 1.4937(1) (zero-density), P ~ 1.3327(1) (FD), to be compared with the 
full-monomer results P ~ 1.40. The state-dependent potential provides now a result which 
is slightly more accurate. 

From these results, it is clear that state-dependent potentials provide only approximations 



to pressure and density that worsen as the fugacity increases. Accurate results can on 



y be 



obtained by using CG multiblob models. If we use the tetramer model of Ref. (42, for 
zR\ = 2.47 we have Z = 2.3457(5) and $ = 0.7978(1), while for zR 3 g = 77.7 we have 
Z = 3.8450(7) and $ = 1.5257(1). Good agreement is observed for the lowest value of the 
fugacity, while small discrepancies occur for the second one. If we use the decamer model^ 
we would obtain for zR g = 77.7 the estimates Z = 3.925(1) and $ = 1.4984(3), which are 
now fully consistent with the full monomer results. 



VI. CONCLUSIONS 

In this paper we have considered structurally derived state-dependent potentials for CG 
models, with the purpose of understanding how useful the y ar e in predicting the thermody- 



namics of complex systems. As already discussed in Ref. |37|, the virial pressure computed 
in the CG model does not reproduce the pressure of the underlying system. But the same 
conclusions apply to essentially all thermodynamic quantities. In the canonical ensemble 
Widom's method does not give the correct underlying-model chemical potential, while in 
the grand canonical ensemble the virial pressure or (N)/V do not reproduce pressure and 
density of the underlying system. We also considered the Ascarelli-Harrison pressure:— it 
is not exact and represents a worse approximation than the virial pressure. It is clear from 
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the discussion that these results are not specific to these quantities. If the pair distribution 
function is used to define the state-dependent potentials, only the quantities that are directly 
related to the pair distribution function, like K(p) or H(z), are identical in the CG and in 
the un derly ing model. All others differ. 

Ref. [37] noted that the correct pressure could be obtained by using the compressibility 
route. However, it should be clear that one is referring to the compressibility route in the 
underlying model, i.e., to Eq. (j3J), not to the compressibility route applied to any of the 
state-dependent CG models. But this means that one is obtaining the pressure from a 
quantity computed in the underlying microscopic model, hence the DD potential plays no 
role. Analogously, in the grand-canonical ensemble, the correct density can be obtained by 
integrating H(z), but again this calculation does not really make use of the model with FD 
potentials. 

It is important to stress that structurally derived state-dependent potentials depend on 
the ensemble one uses, or rather on the thermodynamic variable one uses to identify the 
thermodynamic state of the system. At a given thermodynamic state of the microscopic 
system, DD potentials and FD potentials differ, since they correspond to CG monoatomic 
systems that have the same pair distribution function but different densities and chemical 
potentials. Therefore, DD canonical potentials should not be used to predict phase lines, 
since the coexistence analysis is typically performed in other ensembles. 

Although in this paper we have only discussed DD and FD potentials, the results are 
general and apply to other ensembles, for instance to pressure dependent potentials in the 
isobaric ensemble. However, the conclusions do not necessarily apply to potentials that are 
only temperature dependent. If the thermodynamic state one is considering can be connected 
by a thermodynamic path at fixed temperature to the zero-density state (this is always true 
in the absence of phase transitions or at least in the low-density phase), since temperature 
does not appear explicitly in the integrations leading to the free energies, the free energy is 
correct. Hence, all thermodynamic quantities that can be obtained from the free energies 
are correct. An interesting method for transferring CG potentials between temperatures in 



54. 



the context of the multiscale coarse-graining force-matching method is described in Ref. 
It is important to note that the ensemble dependence does not occur in CG models defined 
by using the multiscale coarse-graining approach based on the force-matching method.— 
Indeed, the distribution of the atomistic configurations is obviously ensemble independent 
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in the infinite- volume limit, hence a 
depend on the ensemble (see Ref. 
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so the matching conditions and the CG model do not 
for a discussion in the isobaric-isothermal ensemble). 
Our results also shed some light on mixed approaches in which one matches the pair 
distribution and some other thermodynamic property. These approaches are conceptually 
correct only if the thermodynamic quantity is appropriately chosen. In the canonical ensem- 
ble, the only correct choice is the isothermal compressibility, i.e. the quantity K we have 
introduced, while in the grand canonical ensemble one should use H(z). It is also possible to 
enforce the value of the pressure, i.e. to work in the isobaric ensemble, but in this case one 
should relax the condition that the microscopic system and the CG model have the same 
density. For each choice, a different state-dependent potential is obtained, as observed, for 



instance, in Ref. [56| in the context of CG models of water. Note that these observations 
can be trivially generalized to multicomponent systems. In this case, one can analogously 
enforce the Kirkwood-Buff 57 equations that relate thermodynamic properties to integrals of 
distribution functions, see Ref. |58| for a recent example. 

The fact that DD potentials are unable to provide the correct results for the thermo- 
dynamic quantities should not, a posteriori, be surprising. After all, even with the use of 
a state-dependent potential, the CG model is still an approximation to the original mi- 
croscopic model. Hence, one should consider this approach as a simple method to obtain 
relatively good approximations for the behavior of complex systems. It is easy to find sit- 
uations, both in theoretical or experimental work, in which state-dependent potentials are 
useful. For instance, one could consider a system of complex molecules for which simula- 
tions are particularly difficult. In order to obtain with reasonable precision the value of the 
thermodynamic quantities in a range of densities, one could perform a run at a density p p , 
measuring all thermodynamic quantities and the density-dependent potential Vd,cg{ y 'i Pp)- 
Then, to compute the pressure for p ^ p p , one could write 

P(p) = P(p p ) + [P(p)-P(p p )}, (75) 

and determine the second term by using the CG model that uses the DD potential computed 
at density p p . Hence, one could estimate 

P(p) = P( Pp ) + [P CG (P, P P ) ~ Pcg(p p , P P )}, (76) 

using the CG model to compute the differences between the state point of interest and 
that for which exact results, i.e. results obtained in the microscopic model, are available. 
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CG models that use zero-density potentials use exactly this strategy, fixing p p = 0. The 
same strategy could also be used in experiments, allowing one to obtain all thermodynamic 
informations in a range of densities, using experimentally determined thermodynamic and 
structural data at p p . 

We have checked this strategy in our polymer example. Assuming that the exact results 
for p and \x are available at $ = 0.8, we wish to compute p and /i at $ = 1.5. Using the 
DD potential computed at $ = 0.8, we obtain Z = 2.44790(4) and /3p exc = 2.8108(3) at 
$ = 0.8 (see Table Q} and Z = 3.8332(4) and /3/i exc = 5.4962(4) at $ = 1.5. Hence, we 
would estimate at $ = 1.5 

Z = 3.74 f3p exc = 5.25, (77) 

to be compared with the FM predictions Z « 3.90 and /3/x exc = 5.38. Estimates (177)) are 
significantly more precise than those reported in Table HI confirming the usefulness of the 
approach. 

Finally, note that we have considered here very simple CG models in which the CG 
molecules are monoatomic. In practical soft-matter applications CG models are usually 
polyatomic. It is quite obvious that the same conclusions apply to these more complex 
models: also in this case structurally derived potentials depend on the ensemble and ther- 
modynamics is only approximately reproduced. 

Acknowledgments 

The authors thank Barbara Capone and Luigi Delle Site for useful comments. CP. is 
supported by the Italian Institute of Technology (IIT) under the SEED project grant number 
259 SIMBEDD Advanced Computational Methods for Biophysics, Drug Design and Energy 
Research. 



' Electronic address: | giuseppe.dadamo@aquila.infn.it 



* Electronic address: andrea.pelissetto@romal.infn.it 



Electronic address: carlo.pierleoni@aquila.infn.it 



23 



1 G. A. Voth, ed., Coarse- Graining of Condensed Phases and Biomolecular Systems (CRC Press, 
Boca Raton, 2009). 

2 R. Feller, Guest Editor, Phys. Chem. Chem. Phys. 11, 1853 (2009). 

3 M. Wilson, Guest Editor, Soft Matter 5, 4341 (2009). 

4 Multiscale Modelling of Soft Matter, Faraday Discussion 144, 1 (2010). 

5 L. Larini, L. Lu, and G. A. Voth, J. Chem. Phys. 132, 164107 (2010). 

6 G. S. Rushbrooke and M. Silbert, Mol. Phys. 12, 505 (1967). 

7 J. S. Rowlinson, Mol. Phys. 12, 513 (1967). 

8 J. A. Barker, D. Henderson, and W. R. Smith, Mol. Phys. 17, 579 (1969). 

9 G. Casanova, R. J. Dulla, D. A. Jonah, J. S. Rowlinson, and G. Saville, Mol. Phys. 18, 589 
(1970). 

M. A. van der Hoef and P. A. Madden, J. Chem. Phys.lll, 1520 (1999). 

1 A. A. Louis, P. G. Bolhuis, J. P. Hansen, and E. J. Meijer, Phys. Rev. Lett. 85, 2522 (2000); 
P. G. Bolhuis, A. A. Louis, J. P. Hansen, and E. J. Meijer, J. Chem. Phys. 114 4296 (2001). 

2 F. Muller-Plathe, Chem. Phys. Chem. 3, 754 (2002). 

3 C. Peter and K. Kremer, Soft Matter 5, 4357 (2009). 

4 S. Izvekov and G. A. Voth, J. Phys. Chem B 109, 2469 (2005); J. Chem. Phys. 123, 134105 
(2005). 

5 W. G. Noid, J.-W. Chu, G. S. Ayton, V. Krishna, S. Izvekov, G. A. Voth, A. Das, and H. C. 
Andersen, J. Chem. Phys. 128, 244114 (2008); W. G. Noid, P. Liu, Y. Wang, J.-W. Chu, G. S. 
Ayton, S. Izvekov, H. C. Andersen, and G. A. Voth, J. Chem. Phys. 128, 244115 (2008). 

6 W. G. Noid, J. W. Chu, G. S. Ayton, and G. A. Voth, J. Phys. Chem. B 111, 4116 (2007); J. 
W. Mullinax and W. G. Noid, Phys. Rev. Lett. 103, 198104 (2009). 

7 S. Merabia and I. Pagonabarraga, J. Chem. Phys. 127, 054903 (2007). 

8 E. C. Allen and G. C. Rutledge, J. Chem. Phys. 130, 034904 (2009). 

9 R. L. Henderson, Phys. Lett. 49A, 197 (1974). 

;0 B. M. Axilrod and E. Teller, J. Chem. Phys. 11, 299 (1943). 
;1 C. N. Likos, Phys. Rep. 348, 267 (2001). 

;2 C. N. Likos, H. Lowen, M. Watzlawek, B. Abbas, O. Jucknischke, J. Allgaier, and D. Richter, 

Phys. Rev. Lett. 80, 4450 (1998). 
:3 M. Watzlawek, C. N. Likos, and H. L6wen, Phys. Rev. Lett. 82, 5289 (1999). 

24 



24 A. Jusufi, J. Dzubiella, C. N. Likos, C. von Ferber, and H. Lowen, J. Phys.: Condens. Matter 
13, 6177 (2001). 

25 J. Dzubiella, C. N. Likos, and H. Lowen, J. Chem. Phys. 116, 9518 (2002). 

26 G. D'Adamo, A. Pelissetto, and C. Pierleoni, J. Chem. Phys. 136, 224905 (2012). 

27 D. Reith, M. Piitz, and F. Muller-Plathe, J. Comput. Chem. 24, 1624 (2003). 

28 A. Lyubartsev and A. Laaksonen, Phys. Rev. E 52, 3730 (1995). 

29 A. Soper, Chem. Phys. 202, 295 (1996). 

30 H. Meyer, O. Biermann, E. Faller, D. Reith, and F. Muller-Plathe, J. Chem. Phys. 113, 6264 
(2000). 

31 R. L. C. Akkermans and W. J. Briels, J. Chem. Phys. 114, 1020 (2001). 

32 J. P. Hansen and I. McDonald, Theory of Simple Liquids, 3rd ed. (Academic Press, Amsterdam, 
2006) 

33 M. S. Shell, J. Chem. Phys. 129, 144108 (2008). 

34 V. Krishna and L. Larini, J. Chem. Phys. 135, 124103 (2011). 

35 V. Riihle, C. Junghans, A. Lukyanov, K. Kremer, and D. Andrienko, J. Chem. Theory Comput. 
5, 3211 (2009). 

36 F. H. Stillinger, H. Sakai, and S. Torquato, J. Chem. Phys. 117, 288 (2002). 

37 A. A. Louis, J. Phys.: Condens. Matter 14, 9187 (2002). 

38 M. E. Johnson, T. Head-Gordon, and A. A. Louis, J. Chem. Phys. 126, 144509 (2007). 

39 P. Ascarelli and R. Harrison, Phys. Rev. Lett. 22, 285 (1969). 

40 B. Widom, J. Chem. Phys. 39, 2802 (1963). 

41 G. D'Adamo, A. Pelissetto, and C. Pierleoni, Soft Matter 8, 5151 (2012). 

42 G. D'Adamo, A. Pelissetto, and C. Pierleoni, J. Chem. Phys. 137, 024901 (2012). 

43 As in our previous work, we determine g(r; p) by using the Domb-Joyce model at w = 0.505838, 



see Ref. 



59l for a precise definition of the model, simulating lattice chains of length 2400. As 



discussed in Ref. [59j this should provide an accurate estimate of the pair distribution function 
in the asymptotic (infinite degree of polymerization), universal limit. 

44 P. G. Bolhuis and A. A. Louis, Macromolecules 35, 1860 (2002). 

45 A. Pelissetto and J.-P. Hansen, J. Chem. Phys. 122, 134904 (2005). 

46 A. Pelissetto, J. Chem. Phys. 129, 044901 (2008). 

47 C. Pierleoni, B. Capone, and J. P. Hansen, J. Chem. Phys. 127, 171102 (2007). 

25 



D. Fritz, V. A. Harmandaris, K. Kremer, and N. F. A. van der Vegt, Macromolecules 42, 7579 
(2009). 

A. J. Clark and M. G. Guenza, J. Chem. Phys. 132, 044902 (2010); A. J. Clark, J. McCarty, I. 
Y. Lyubimov, and M. G. Guenza, Phys. Rev. Lett. 109, 168301 (2012). 
T. Vettorel, G. Besold, and K. Kremer, Soft Matter 6, 2282 (2010). 
T. Morita, Prog. Theor. Phys. 23, 839 (1960). 
P. Attard, J. Chem. Phys. 94, 2370 (1991). 

In the full-monomer model we have (3PRg = 0.45 and (3PRg = 1.40 for = 0.80 and 1.50. For 

comparison, we report the pressure-dependent potentials for these two values of the pressure. 

Using the HNC approximation we can also compute the thermodynamic variables using these 

state-dependent interactions. We obtain $ = 0.784, Z = 2.397, /3/i cxc = 2.72 for f3PR 3 g = 0.45 

and $ = 1.484, Z = 3.942, /3^ exc = 5.71 for ftPR^ = 1.40. 

V. Krishna, W. G. Noid, and G. A. Voth, J. Chem. Phys. 131, 024103 (2009). 

A. Das and H. C. Andersen, J. Chem Phys. 132, 164106 (2010). 

H. Wang, C. Junghans, and K. Kremer, Eur. Phys. J. E 28, 221 (2009). 

J. G. Kirkwood and F. P. Buff, J. Chem. Phys. 19, 774 (1951). 

P. Ganguly, D. Mukherji, C. Junghans, and N. F. A. van der Vegt, J. Chem. Theory Comput. 
8, 1802 (2012). 

S. Caracciolo, B. M. Mognetti, and A. Pelissetto, J. Chem. Phys. 125, 094903 (2006). 



26 



